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ABSTRACT 



Aims. We present the current status of ongoing searches for molecular hydrogen in high-redshift (1.8 < z a bs S 4.2) Damped Lyman-o- 
systems (DLAs) capitalising on observations performed with the ESO Very Large Telescope (VLT) Ultraviolet and Visual Echelle 
Spectrograph (UVES). 

Methods. We identify 77 DLAs/strong sub-DLAs, with logiV(Hi) > 20 and z abs > 1.8, which have data that include redshifted H 2 
Lyman and/or Werner-band absorption lines. This sample of H i, H2 and metal line measurements, performed in an homogeneous 
manner, is more than twice as large as our previous sample (Ledoux et al. 2003) considering every system in which searches for H 2 
could be completed so far, including all non-detections. 

Results. H2 is detected in thirteen of the systems, which have molecular fractions of values between / = 5 x 10~ 7 and / = 0.1, 
where / = 2AT(H2)/(2Af(H 2 ) + A'(Hl)). Upper limits are measured for the remaining 64 systems with detection limits of typically 
log./V(H 2 ) ~ 14.3, corresponding to log/ < -5. We find that about 35% of the DLAs with metallicities relative to solar [X/H] > —1.3 
(i.e., l/20 th solar), with X = Zn, S or Si, have molecular fractions log/ > -4.5, while H 2 is detected - regardless of the molecular 
fraction - in ~ 50% of them. In contrast, only about 4% of the [X/H] < -1.3 DLAs have log / > -4.5. We show that the presence of 
H2 does not strongly depend on the total neutral hydrogen column density, although the probability of finding log / > -4.5 is higher 
for logA'(Hl) > 20.8 than below this limit (19% and 7% respectively). The overall H 2 detection rate in logJV(Hl) > 20 DLAs is 
found to be about 16% (10% considering only log / > -4.5 detections) after correction for a slight bias towards large A'(Hl). There 
is a strong preference for H 2 -bearing DLAs to have significant depletion factors, [X/Fe] > 0.4. In addition, all FL-bearing DLAs have 
column densities of iron into dust grains larger than log Ntpe)^ ~ 14.7, and about 40% of the DLAs above this limit have detected 
H 2 lines with log/ > -4.5. This demonstrates the importance of dust in governing the detectability of H 2 in DLAs. Our extended 
sample supports neither the redshift evolution of the detection fraction of H 2 -bearing DLAs nor that of the molecular fraction in 
systems with H 2 detections over the redshift range 1.8 < z a bs 5 3. 

Key words, cosmology: observations - quasars: absorption lines - galaxies: ISM - ISM: molecules 



1. Introduction 

Damped Lyman-a systems (DLAs) we r e discovered in the 
seventies (e.g . . lLowrance et al. Il972t iBeaver etafj Il972t 
ICarswell et alj Il975t IWright et al 



of present-day galaxies (see, e.g., Wolfe & Prochaskal [2000; 



e.g.. iwo 

lHaehneltetai]|2000l IWolfe et alj|2005t) . Our understanding of 
DLAs is mainly b ased on the study of low- ionisation metal 
absorptions (e.g., Prochaska & Wolfe|_ 20021) but also high- 
ionisation specie s dLu et alj 119961 ; IWolfe & Prochaskal [2000; 
iFox et al. Il2007allbh and, i n a few cases, molecular absorptions 
(e.g., Ledoux et al.ll2003h . The latter are not conspicuous how- 
ever in contrast to what is seen in the Galaxy and, for a long time, 
only the D LA towards Q 0528-2505 was kn own to contain H2 
molecules dLevshakov & Varshalovichll985l) . H2 -bearing DLAs 
are nevertheless crucial to understand the nature of DLAs be- 
cause molecular hydrogen is an importa nt species to derive the 
physical conditions in the gas (see, e.g.. | Tumlinson et al.ll200: 
i Reimers et alj l2003t iHirashita & FerraraT 20051; ISrianand et 
l2005HCui et al.ll2005HNoterdaeme et al.ll2007bl) 

The first systematic search for molecular hydrogen in high- 
redshift (z a b s > 1 .8) DLAs was carried out using the Ultraviolet 
and Visual Echell e Spectrograph (UV ES) at the Very Large 
Telescope (VLT) dLedoux et al.l 12003). It consisted of a sam- 
ple of 33 DLAs with H2 detected in eight of them. Molecular 



1979) and identified after- 



wards as redshifted damping absorptio ns from large c olumn den- 
sities of neutral atomic hydrogen (Smit h et"aT1l 1979h . 

Numerous DLAs, with JV(Hi) > 2 x 10 20 atoms ctrT 2 , 
have been discov ered through large dedicated surveys (e.g., 
IWolfe etal.lll986h and more recently thanks to the huge num- 
ber of quasar spectra availab le from the Sloan Digital Sky 
Survey dProchaska et al.l 120051) . Because DLAs contain most 
of the ne utral hydrogen available for star form ation in the 
Universe (Wolfe et al. 1986; lEanzetta et all 1 199 lh and are as- 
sociated with numerous metal absorption lines, they probably 
arise in the interstellar medium of protogalaxies, progenitors 



* Based on data gathered at the European Southern Observatory 
(ESO) using the Ultraviolet and Visual Echelle Spectrograph (UVES) 
installed at the Very Large Telescope (VLT), Unit-2 (Kueyen), on Cerro 
Paranal, Chile. 
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fractions were found to lie in the range -3.5 < log/ < -1 
with / = 2JV(H 2 )/(2jV(H 2 ) + N(H i)). Upper limits of typically 
Af(H 2 ) ~ 2 x 10 14 cm' 2 (corresponding to log/ < -5) were 
measured in the other systems. More recently, we noted a cor- 
relation between the presence of molecular hydrogen a nd the 
metallicity of high-redshift DLAs dPetitjean et al.ll2006l) . High 
molecular fractions (log / > -4) were found in about 40% of 
the high-metallicity DLAs ([X/H] > l/20 th solar) whilst only 
~ 5% of the [X/H] < -1.3 DLAs have log/ > -4. Other 
papers by our group focused on specific detections. We pre- 
sented the analysis of three systems with low molecular frac- 
tions, i.e., log/ < -4, o ne of them having a low metallicity 
dNoterdaeme et al.l l2007al) . and the H 2 -bearing DLA with, to 
date, the highest reds hift, at z abs = 4.224 towards Q 1441+2737 
dLedoux et al.l2006bh . We note that a possible detection of H 2 in 
a DLA towards a Gamma-ray Burst ( GRB) afterglow has been 
reported recently dFvnbo et al.l 1200 6). However, the origin of 
DLAs at the GRB host-galaxy redshift is very li kely to be differ- 
ent from those observed in Q SO spectra (e.g., iJakobsson et al.l 
120061; iProchaska et alj|2007ab . 

We present here the whole sample of UVES high-redshift 
QSO-DLAs for which the wavelength range where H 2 lines are 
redshifted is covered by the available spectra. T his sample is 
more than twice as large as in our previous study dLedoux et al.l 
120031) . We present the observations and the UVES DLA sam- 
ple in Sect. 2 and provide comments on individual absorbers in 
Sect. 3. We discuss the overall population in Sect. 4 and results 
in Sects. |5]to[8] We conclude in Sect. [9] 

2. Observations and sample 

All the quasars in our sample were observed with the Ultraviolet 
and Visual Echelle Spectrograph (UVES, [Dekkeret all [2000) 
mounted on the ESO VLT-UT2 (Kueyen) 8.2 m telescope on 
Cerro Paranal, Chile. We have used our UVES database to build 
up a sample of 77 DLAs/strong sub-DLAs along 65 lines of 
sight, hereafter called sample S. The systems were selected as 
having jV(Hi) > 10 20 cirT 2 , and redshifts z abs > 1.8 so that at 
least part of the wavelength range into which H 2 absorption lines 
are expected to be redshifted, is covered. Systems for which the 
flux at the wavelengths of all redshifted H 2 lines is zero as a 
consequence of the presence of an intervening Lyman-limit sys- 
tem located at a higher redshift were rejected. Systems for which 
the Ly-ff forest is so crowded that no meaningful upper limit on 
Af(H 2 ) could be derived were also excluded from the final sam- 
ple. 

Sample S comprises a total of 68 DLAs (logiV(Hl) > 20.3) 
and nine strong sub-DLAs (20 < log N(H i) < 20.3) according to 
the co mmon definition of DLAs (log jV(H i) > 20.3; Wol fe et all 
119861) . Our lower limit on N(Hi) for sub-DLAs is only slightly 
lower than the classical definition of DLAs ensuring that these 
absorbers are mostl y neutral and share the same physical nature 
as classical DLAs dViegaslll995l) . 

Most of the systems in sample S (53 out of 77) come from 
the sample of iLedoux et al.l (l2006al) . This sample is mainly 
drawn f rom the follow-u p of the Large Bright QSO Survey 
(LBQS, IWolfeetal.1 1995b and has been observed between 2000 
and 2004 in the course of our systemati c search for molec- 
ular hydrogen at z abs > L8 (see also iPetitiean et al.1 120001; 
ILedoux et al] 120031) . This comprises 46 bona-fide DLAs and 
seven strong sub-DLAs. 

In addition to this sample, we are considering UVES data 
for 13 absorbers (11 DLAs a nd two strong sub-D LAs) from the 
Hamburg-ESO DLA survey dSmette et al.l 120051) . seven DLAs 



from the CORALS survey (lAkerman et al.l 120051) . and four 
DLAs mainly from our own observing runs. Among the latter 
systems, the Zabs = 3.692 and3.774DLAstowardsQ0131+0345 
and the z abs = 2.659 DLA towards Q 0642-5038 were ob- 
served in visitor mode on September 17-20 2004, under Prog. 
ID 073.A-0071 (PI: Ledoux). On the other hand, thez abs = 4.203 
DLA towards Q 095 1-0450 was observed in service mode on 
January 26-27 and February 19-21, 2004, under Prog. ID 072.A- 
0558 (PI: Vladilo). 

We have reduced all the data including those retrieved from 
the ESO archive in an homo geneous manner using the UVES 
pipeline (Ball ester et al.l 120001) . which is available as a dedi- 
cated package of the ESO MIDAS data reduction system. The 
main characteristics of the pipeline are to perform a robust inter- 
order background subtraction for master flat-fields and science 
frames and an optimal extraction of the object signal subtract- 
ing the sky spectrum and rejecting cosmic rays simultaneously. 
The wavelength scale of the reduced spectra was converted to 
vacuum-heliocentric values. Each spectrum, corresponding to 
different instrument settings, was rebinned to a constant wave- 
length step. No further rebinning was performed during subse- 
quent data analysis. Individual exposures were then weighted, 
scaled and combined altogether. 

Total neutral hydrogen column densities have been mea- 
sured from Voigt-p rofile fitting of the d amped Lyman-o' and/or 
Lyman-/? lines (see [Ledoux et al]|2006al Smette et al., in prep). 
Most of the metalli city and depletion measurements are from 
ILedoux et al.l d2006al) except when more recently re-measured 
(see details in footnote of Table [TJ. Apart from the DLAs to- 
wards Q 01 3 1+ 0345 and Q 095 1-0450 (for which abundances 
were taken from IProchaska et all (|2007b) because UVES obser- 
vations did not allow a sufficient number of metal lines to be 
covered), we (re)determined all metal column densities by the 
homogeneous fitting of Voigt-profiles to non-saturated absorp- 
tion lines. This includes the systems from the Hamburg-ESO 
and CORALS DLA surveys. Metallicities were determined us- 
ing the zinc abundance (when Znn is detected), or from the 
abundances of sulphur or silicon. No ionisation correction has 
been applied. Zinc and sulphur are known to be little depleted 
into dust grains. Silicon, in turn, is probably mildly depleted and 
has been used only in cases where Zn and S are not detected. The 
results are summarised in TableQ] Abundances are given relative 
to solar, i.e., [X/H ] = log N(X)/N (H) - log (X/H) . Solar abun- 
dances as listed inlMortonl d2003l) . based on meteoritic data from 
Grevesse & Sauvall (120021) . were adopted. 

3. Comments on individual systems 

Q 0000-2619, z a bs = 3.390: iLevshakov et alJ d2000l 1200 lh re- 
ported a tentative detection of H 2 in this system. However, they 
detected only two weak absorption features in the Lyman-or for- 
est, identified as H 2 W2-0Q1 and L4-0R1. The probability that 
these features are actually due to intervening Lyman-o' absorbers 
is high and we consider the derived 7=1 column density as an 
upper limit only. In any case, this system has a very low molec- 
ular fraction of log / < -6.87. 

0013-0029^^ = 1.973: The detection of H 2 in this sys- 
tem h a s been reported for the first time by iGe & Bechtoldl 
dl997h . iPetitjean efall (120021) showed that this system is actu- 
ally the blend of a DLA (logJV(Hi) = 20.83) and a sub-DLA 
(logA^(Hi) < 19.43) separated by ~500kms _1 .H 2 is present in 
both systems, with four components detected up to J = 5. 
Q0027-1836, z abs = 2.402: This system has the lowest metal- 
licity ([X/H] = -1.63) of systems in which H 2 molecules have 
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been detected. H2 is detected in a single component up to rota- 
tional level 7 — 5 and possibly 6. Thanks to the high data qual- 
ity, we detected an increase in the Doppler parameter b, from 
low to high rotational lev els, for the first time at high redshift 
dNoterdaeme et al.ll2007al) . 

Q 0347-3819, z abs = 3. 025: The detectio n of H 2 in this sys- 
tem was first reported by Levsha kov et al.l d2002l) using UVES 
commissioning data. Subsequently iLedoux et alj d2003l) anal- 
ysed higher quality spectra of the system. Weak molecular ab- 
sorption is observed in a single component up to 7 = 4. 
Q 0405-441 8, = 2.595: Molecular h ydrogen is detected 
up to 7 - 3 in a s ingle component (see ILedoux et al.l 120031: 
ISrianand et al.ll2005b . There is some indication that the relative 
strengths of the 7 — 2 and 3 lines require b S 1.5 kms -1 , while 
the 7 — and 1 lines are consistent with b — 1 .3 km s _1 . 
Q 0528-2505, z abs = 2 811: This is the first DLA syste m in 
which H2 was detected (Levshakov & Vars halovichlll985j). The 
detect i on of H2 was readdre s sed in mo re detail by [F oltz et al. 
(1 19881) . ISrianand & Petitjeanl dl998h and lSrianand et al.l (12005). 
The singl e H2 component seen in a low resolution CASPEC 
spectrum ISrianand & Pe titiean 1998) i s resolved in to two com- 
ponents in a UVES spectrum dSrianand et al.ll2005l) . Rotational 
levels up to 7 = 5 are detected. 

Q0551-3638, z a bs = 1.962: This system has the highest metal- 
licity in the UVES DLA sample with [Zn/H] = -0.35. H 2 is 
clearly det ected in two well-re solved features separated by about 
55 kms 1 dLedoux et al.ll2002h . The first feature is narrow and 
weak, but detected up to 7 = 3. The second feature is broader 
and stronger and has been fitted using two components at pre- 
cisely the same redshift as the detected C 1 lines. 
Q 0642-5038, z a bs = 2.659: FL- in this system is detected in a 
single strong component. This is a new detection and details will 
be presented in a forthcoming paper (Ledoux et al ., in prep.). 
Q0841 + 1256, z abs = 2.375: lPetitjean et al.l d2000l) detected two 
weak absorption features corresponding to H2 L4-0 R0 and L2- 
R0 at the same redshift. As for the DLA towards Q 0000-26 19, 
we consider the tentative measurement of the 7 = column den- 
sity as an upper limit only. 

Q 1 232+ 08 1 5, z a bs = 2.33 8: A single H2 component is seen i n 
rotational levels 7 = to 5 dGe et al.l2001b ISrianand et al.l2 005). 
This is the only DLA in which de uterated molecular hydr ogen 
(HD) has been detected until now dVarshalovich et alJ UOQl). 
Q 1441+2737, z abs = 4.224: Because of the high redshift of this 
system (actually the highest amongst known ^-bearing DLAs), 
the Lyman-or forest is quite dense in the wavelength region in 
which H2 lines are redshifted. However, H2 could be detected 
and studied in a large number of transitions thanks to the high 
spectral resolution of the data. Four rotation al levels are detecte d 
in no less than three velocity components dLedoux et a l. 2006b). 
Q 1444+0126, z a b s = 2.087: Molecular hydrogen is found in 
two components of this sub-DLA in 7 < 3 rotational levels 
dLedoux et all2003l) . 

Q 1451 + 1223, Zabs = 2.469: Because of the Lyman-break from 
the DLA at z a bs = 3.171, the expected positions of only two lines 
(H2 L0-0R0 and L0-0R1) fall in a region of non-zero flux. The 
high upper limit we derive is not strongly constraining and the 
molecular fraction is probably well below log / = -4.5. 
Q2318-1107,z abs = 1.989: Weak H 2 absorption features from 
7 = to 2 are una mbiguously detected in a single component 
dNoterdaeme et alJl2007al) . 

Q 2343+ 1232, z abs = 2.431: Thanks to the high quality of the 
spectrum, H2 absorption lines (from 7 — and 1) are unam- 
biguously detected although they are extremely weak (A^(H2) ~ 
5 x 10 13 cm -2 ). We measure the lowest molecular fraction yet 



measured for a DLA, l og/ = -6.41 dPetitiean et al.l l2006t 
iNoterdaeme et al.ll2007al) . 

Q 2348-0108, z a bs = 2.426: No less than seven H 2 com- 
ponents are detected in this system (Petiti ean et al.l [2006; 



INoterdaeme et al.ll2007bl) . The H2 profile is complex with three 
strong and four weak components spread over about 250 km s _1 , 
making H2 lines from different transitions overlap. However, the 
large number of observed transitions and the good data quality 
allow an accurate measurement of column density for different 
rotational levels. 



4. Overall population 

In this Section, we use the results summarised in Table [TJ to de- 
rive the properties of the H2 gas in the global DLA population. 
We list in this table, the name of the quasar in column 1, the 
emission and absorption redshifts in columns 2 and 3, the total 
H 1 column density in column 4, the metallicity relative to solar 
[X/H] and the depletion factor [X/Fe] in columns 5 and 6, with 
X given in column 7, t he velocity spread of the low-ionisation 
metal line profiles (see ILedoux et al.ll2006al) in column 8, the 
total molecular hydrogen column densities in the 7 = and 1 
rotational levels in columns 9 and 10, the mean molecular frac- 
tion in column 11, and the references for the measurements in 
column 12. Upper limits on the molecular fraction were calcu- 
lated to be the sum of the upper limits on the first two rotational 
levels (7 = and 1) and are given at the 3 <x significance level. 
The thirteen firm detections of H2 constitute what we will here- 
after call sub-sample Shi- Ten of them have molecular fractions 
higher than all upper limits, i.e., log / > -4.5, and will consti- 
tute sub-sample Shf (standing for high molecular fractions). We 
note that obviously Shf c Shi c S. 

The upper panel of Fig. [TJ gives the distribution of A^(H 1) 
from sample S, sub-samples Shi and Shf to be compared to 
the reference distribution derived from the SDSS-DR5 DLA 
samp le scaled to the num ber of systems in the UVES sam- 
ple dProchaska et al.ll2005l) . It is apparent that the distribution 
from the UVES sample is slightly biased compared to the scaled 
SDSS distributiorQ in favour of large N(U 1) DLAs (i.e., with 
logAf(Hi) > 20.8). We can however correct for this bias. We 
divide each DLA sample into two halves, above and below 
logAf(Hi) = 20.8. We scale the number of systems in each bin 
of the UVES distribution by the ratio of the number of systems 
in the (scaled) SDSS and UVES samples, above (by a ratio of 
~ 0.5) and below (by ~ 1.4) log N(H 1) = 20.8. The result of this 
scaling is shown in the bottom panel of Fig. [TJ In the inset, the 
cumulative distributions are shown in order to compare the two 
populations. It is apparent that the scaling has corrected for the 
above-mentioned bias. 

The overall H2 detection rate in log A^(H 1) > 20 systems is 
found to be about 16% (10% considering only log/ > -4.5 de- 
tections) after correction for the above-mentioned bias. There is 
no H2 detection for log A^(H 1) < 20.2. These numbers, however, 
may well be affected by small number statistics. H2 is in fact 
detected in a sub-DLA component (log A^(H 1) ~ 19.4) of the 
Zabs = 1-973 system towards Q 0013-0029 (see iPetitiean et al.l 
|2002|) . The corresponding molecular fraction is larger than 10~ 4 . 
In any case, additional observations are required to investigate 
the molecular content of sub-DLAs. Apart, possibly, from the 
first bin (log N(H 1) < 20.2), there is no strong dependence of 



1 Note that the SDSS A'(Hi) distribution itself may be biased, how- 
ever it is currently the best available representation of the DLA popula- 
tion. 
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the presence of molecules on the neutral hydrogen column den- 
sity. Large molecular fractions are seen over the range 20.2 < 
log 7V(H i) < 21.8. 
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Fig. 1. Top panel: Neutral hydrogen column density distribu- 
tions of DLAs in the overall UVES sample S (solid), the sub- 
sample of H2-bearing systems Shi (grey), and that of the sys- 
tems with log / > -4.5 (.Shf; d ark grey). The distribu tion from 
the SDSS-DR5 sample (dotted; iProchaska et al.ll2005l) is repre- 
sented with a different scaling (right axis) so that the area of both 
histograms are the same. The inset shows the cumulative distri- 
butions from sample S (solid line), sub-sample Shf (dashed), 
and the SDSS-DR5 sample (dotted). Bottom panel: The above 
log N(H i) distributions are corrected for the bias in favour of 
large N(H i) column densities (see text) by scaling the numbers 
separately for 20 < logJV(Hi) < 20.8 and log N(Hi) > 20.8 to 
match the SDSS-DR5 distribution. 



The double-sided Kolmogorov-Smirnov test yields a proba- 
bility Pks = 0.57 that sample S and sub-sample Shf are drawn 
from the same parent population. 

However, molecular hydrogen with log / > -4.5 is detected 
in respectively 9% and 19% of the systems in the two sub- 
samples: 20.2 < logJV(Hi) < 20.8 (3/35) and logJV(Hi) > 20.8 
(7/36), which implies that there is probably a tendency for H2 to 
be seen more often in large column density systems. We note that 
the fraction is 7% (3/41) for the range 20.0 < log/V(H 1) < 20.8. 
We also note that we are considering here the total neutral hy- 
drogen column densities. The actual H i-clouds corresponding to 
the H2 absorptions probably have smaller N(H 1). This is in fact 
observed for the DLA atz abs = 1.973 towards Q 0013-0029, for 
which partial deblending of the Ly-a absorption is possible. 



The fact that the sample is biased in favour of large Hi 
column density DLAs could be a problem for the analysis of 
the overall DLA population if any correlation exists between 
log/Y(Hi) and the metallicity ([X/H]) and/or the depletion of 
metals onto dust grains ([X/Fe]). This is not the case how- 
ever: Figure [2] shows the neutral hydrogen column density as a 
function of metallicity (top panel) and depletion factor (bottom 
panel). It is apparent that there is no correlation between [X/H] 
or [X/Fe] and logTY(Hi). The Kendall's rank correlation coeffi- 
cient is as small as 0.02 (resp. -0.01) for [X/H] (resp. [X/Fe]) vs. 
log N(H 1). The lack of systems with both a high metallicity and 
large N(H 1) is a feature common to all DLA samples. It can be 
explained either by a lack of la rge H 1 colum n densities because 
of the transition of H 1 into H2 dSchavell200ll) and/or by the fact 
that any quasar located behind such an absorber would remain 
undetected because of the ind uced extinction (e.g., iBoisse et al.l 
1 1 9981 : IvTadilo & Perouxl2005l) . Such systems could also be asso- 
ciated with regions with low projected cross-sections. Therefore, 
although the lack of high-metallicity systems with large N(H 1), 
together with the bias towards large N(H 1) in our sample, could 
introduce a small bias towards low metallicities, the absence of 
the above-described correlations implies that any bias will have 
little influence on the properties based on the metallicities and 
depletion factors discussed in Sect. [6] 
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Fig. 2. Logarithm of the total neutral hydrogen column density 
in the UVES sample versus metallicity (top panel) or depletion 
factor (bottom panel). Filled squares indicate systems in which 
H2 is detected. The lack of systems with both a high metallicity 
and a large H 1 column density is seen in all DLA samples. There 
is no correlation between N(H 1) and [X/H] or [X/Fe]. 
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In Fig. [3] we compare the metallicity distribution of DLAs 
in the UVES sample to that of DLA s in the Keck sample 
(HIRES+ESI, iProchaska etail l2007bl) . Although a difference 
can be seen in the range -1.6 < [X/H] < -1, partly due to the 
bias towards large N(H i), no systematic bias is observed for the 
UVES sample, apart from any also affecting the Keck sample, 
as shown by the cumulative distributions (inset). We note that 
the observed distributions are similar in shape to the distribution 
derived using simulations bv lHou et alj (120051 their fig. 3). The 
metallicity distribution for H2-bearing DLAs is also plotted. It is 
apparent from the figure that the distribution of Lb-bearing sys- 
tems is strongly skewed towards high metallicities (see Sect.|6}. 
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Fig. 3. The metallicity distribution of DLAs in the overall UVES 
sample (S; solid) is compared to that from the Keck sample 
(dotted; IProchaska et al.ll2007bl) . adequately scaled (right axis) 
so that the area of both histograms are the same. The distribu- 
tions are similar. Distributions from sub-samples Sm (grey) and 
Shf (dark grey) are also shown. It is clear from this and the 
Kolmogorov-Smirnov test probability (inset) that the distribu- 
tions from EL-bearing systems are skewed towards high metal- 
licities. 
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Fig. 4. Molecular fraction versus total hydrogen (H 1+H2) col- 
umn density for sample S (the filled squares correspond to H2 
detections, i.e., sub-sam ple ^SmX lines of sig ht through the SMC 



(blue dots, top panel; Tumli nson et al.ll2002h. the LMC (red tri- 
angles, middle panel; iTumlinson et al . 2002) an d the G alactic 



disk (green diamonds, bottom panel; ISavage etail 19771) . Upper 
limits are marked by unfilled symbols with downward arrows. 



5. Molecular fraction 

As can be seen in Fig.|4l there is no sharp transition in the molec- 
ular fraction of DLAs at any total hydrogen (H 1+H2) column 
density A^(H), in con trast to what is obse rved at log A^(H) = 20.7 
in the Galactic disk (Savage et al.lll977 , bottom panel of Fig. |4|l 
or at log Af(H) = 20. 4 along high-latitude Galactic lines of sight 
dGillmon et al.ll2.0d6b . We note also that no sharp transition from 
H 1 to H2 is observed in the Magellanic clouds (upper and mid- 
dle panels of Fig. although SMC lines of sight probe larger 
gas column densities. Such transitions are expected to occur at 
logjV(H) > 21.3 in the LMC and logJV(H) > 22 in the SMC 
( ITumlinson et aljr2002f) . The molecular fractions of DLAs are 
similar to values measured for the Magellanic clouds, but are 
lower than measurements for the Galactic disk. The absence 
of a sharp H 1/H2 transition does not contradict the predictions 
by ISchavd J200l) bee ause the co vering factor of large A^EL. ) 
systems is expected to be small dZwaan & Prochaskal 12006). 
Moreover, the induced extinction could be large (see below) and 
such systems may have so far been missed. We should compile 
a representative sample of strong DLAs (logA^(Hi) > 21.5) in 
order to address this point. 



6. The importance of dust 

Figure [5] is an extension of fig. 3 from Petitiean et alj d2006l) to 
the whole UVES sample, with about 2.5 times more systems 
with [X/H] < -1.3. It can be seen from this figure that about 
35% of the [X/H] > -1.3 (l/20 th solar) systems have molecular 
fractions log/ > -4.5, while H2 is detected - regardless of the 
molecular fraction - in ~ 50% of them. In contrast, only about 
4% of the [X/H] < -1.3 DLAs have log/ > -4.5. We remind 
the reader that we use the limit log / > -4.5 to define systems 
with high molecular fraction because all of our upper limits are 
below this value. The lowest metallicity at which H2 has been 
detect ed to date in a DLA is [Zn/H] = -1.63 dNoterdaeme et alj 
2007a]). The corresponding system, at z aDS = 2.402 towards 
Q 0027-1836, has a molecular fraction of log / = -4.15. 

This trend in metallicity supports the idea that the presence 
of dust is an important ingredient in the formation of H2. The 
correlation bet ween the depletion o f metals into dust grains and 
the metallicity dLedoux et alj|2003l) has now been confirmed us- 
ing a sample more than twice the size, and at higher confidence 
level (5.3 cr compared to Act) (see Fig. [SJ. High-metallicity 
DLAs are usually more dusty and, as a consequence, the for- 
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mation rate of H2 onto dust grains is enhanced, while the photo- 
destruction rate is lowered by dust- and self-shielding. 
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Fig. 5. Logarithm of the molecular fraction, / = 
2jV(H 2 )/(2jV(H 2 ) + N(Hi), versus metallicity, [X/H] = 
log N(X)/N(H) - log(X/H) with X=Zn, S or Si. Filled 
squares indicate systems in which H 2 is detected. The horizontal 
(resp. vertical) dashed line indicates the separation between 
what we somewhat arbitrarily call high and low molecular 
fractions, log/ = -4.5 (resp. high and low metallicities, 
[X/H] = -1.3). 
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Fig. 6. A correlation (5.3 <x confidence level) exists in the overall 
UVES sample between dust-depletion factor, [X/Fe], and metal- 
licity, [X/H], with X either Zn, S or Si. 



Figure [7] shows the distribution of depletion factors (defined 
as [X/Fe], with X=Zn, S or Si) for the different samples of sys- 
tems defined in this paper. It is clear from this histogram that H 2 
is found in DLAs with the highest depletion factors. The prob- 
ability that sample S and sub-sample Sm are drawn from the 
same parent population is very small, i.e., Pks ~ 10~ 3 . 

To test the influence of the presence of dust on the 
molecular hydrogen column density, it is worth comparing 
Af(H 2 ) to the column density of iron int o dust, Af(Fe ) d us t = 
(1 - 10- [X/Fel )MX)(Fe/X) dla (see Fig. [Hi. Ivladilo et alj (2006) 
showed that absorbers associated with large column densities of 
iron in dust, generally induce more extinction (i.e. have larger 
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Fig. 7. Distribution of depletion factors from the overall UVES 
sample (S; solid), sub-samples Sm (grey) and Shf (dark grey). 
It is clear from these histograms that the distribution of depletion 
factors for H 2 -bearing systems is different from that of the over- 
all sample. This shows clearly that H 2 -bearing DLAs are more 
dusty than the rest of the DLAs. The probability that the two 
samples are drawn from the same parent population is indeed 
very small: P ks (S/<S H f) ~ 10~ 3 . 



values of Ay). This correlation is consistent with that observe d 
along Galactic interstellar lines of sight (e.g., ISnow et al.l2 002). 
We note that the above expression of /V(Fe)d us t assumes that X 
is a n on-refractory e lement and that the intrinsic [X/Fe] ratio is 
solar. IVladilol (120021) introduced corrections for zinc depletion 
and possible non -solar intrinsi c abundances (i.e., (Fe/X)di a + 
(Fe/X) Q ; see also IVladilo et al.l l2006) that we do not apply here 
because the corrections are small, and not necessary for our 
purposes. This effect is most important for Silicon that can be 
depleted into dust grains, and the depletion factor is underesti- 
mated when using [Si/Fe] . The effect is however probably less 
than 0.3 dex (see lPetitiean et al. 2002) and should have no con- 
sequence on the characteristics of the whole population. When 
[X/Fe] is negative, which happens for three DLAs, A r (Fe)d lls t can- 
not be computed directly. We thus estimated a 3 <x upper limit on 
Af(Fe)dust by considering the upper bound provided by the 3 cr er- 
ror on [X/Fe] . 

It can be seen on Fig. [8] that all H 2 -bearing DLAs have 
^V(Fe)dust > 5 x 10 14 cm -2 , supporting the idea that dust is 
an important ingredient in the formation of H 2 . We note how- 
ever that this colu mn density corresp onds to a small extinc- 
tion (log Ay 1.5: IVladilo et alJ l2006) possibly explaining the 

absence of detectable H 2 in three DLAs at A^(Fe)d ust > 3 x 
10 15 cm" 2 . 

An additional reason why H 2 is not detected in these systems 
could be related to the particle density in the neutral gas being 
too small. The equilibrium between formation and destruction 
of H 2 molecules can be written as Rmi(H) = ^diss"(H 2 ), where 
^diss is the photo-dissociation rate, R is the formation rate, n 
the particle density, n(H) ^ n the proton density and n(H 2 ) the 
H 2 density. By multiplying this expression by the longitudinal 
size of the cloud, one can see that the H 2 column density de- 
pends linearly on the particle density. Large dust column densi- 
ties may occasionally result from large and diffuse clouds where 
the particle density is insufficiently large for H 2 to form effi- 
ciently enough to be detectable. The column densities of iron 
in dust are in fact large (A^(Fe)d us t > 3 x 10 15 cur 2 ) in three 
DLAs without detectable H 2 mainly because of the large neu- 
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tral hydrogen column densities (i.e. log N(VL i) = 21.70 towards 
Q0458-0203; logTV(Hi) = 21.80 towards Q 1157+0128; and 
log N(H I) = 21 .40 towards Q 1209+0919). 

The fact that H2 is detected in some DLAs at low extinction, 
could imply that the weak (due to high depletion) iron compo- 
nent, associated with the FL-bearing c omponent, is hidden in the 
overall metal-line profile (see, e.g., Ledo ux et al.l 12002). This 
may bias the measurements of metallicities and depletion fac- 
tors, and hence the measurement of column density of iron into 
dust grains. In the future, we propose that this issue is more 
deeply investigated. 




14 

log N(Fe iu J 

Fig. 8. H2 column density, NQI2), versus column density of 
iron into dust grains, -/V(Fe)d u . st . The horizontal dotted line at 
log MEL:) - 18.2 (resp. 14.3) shows the median log NQI2) (resp. 
the median upper limit on logA^FL)) for the FL-bearing DLAs 
from sub-sample Sm (resp. the H2-undetected systems from 
sub-sample S - Sm)- There is a clear gap of about four orders 
of magnitude between the two values. The horizontal bar in the 
upper-left corner of the plot shows the typical uncertainty on the 
measured log iV(Fe)d us t- 



at Zabs — 2.43 towards Q 2348-0108 where no less than seven 
H2 co mponents are spread over 250 kms 1 dNoterdaeme et al.l 
l2007bh . and the z abs = 1.973 system towards Q0013-0029 
where two p airs of H2 componen ts are separated by more than 
500 kms" 1 dPetitiean et alj|2002l) . The latter could be consid- 
ered as the blend of tw o DLA systems (see Sect. [3). In addi- 
tion, iFox et al] d2007al) interpreted DLAs where the C iv total 
line width exceeds the escape velocity to be associated with out- 
flowing winds. Using this criterion, there is no indication that H2 
is found preferably in DLAs with outflows. Indeed, considering 
the 54 DLA syste ms common to the present sample and that of 
IFox et all d2007al) . we can estimate that FL is detected in ~16% 
(3/19) of the DLAs with outflows while this percentage is ~22% 
(12/54) for the entire DLA population. 
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Fig. 9. Distribution of th e velocity spread, A V, of low-ionisation 
metal line profiles (see iLedoux et al.ll2006al for details on the 
technique) from the overall LIVES sample (S, 77 systems, solid 
histogram), systems with detected molecular hydrogen (Sm, 13 
systems, grey), and systems with log/ > -4.5 (Shf, ten sys- 
tems, dark grey). It is apparent from this and the Kolmogorov- 
Smirnov test probability (inset) that the likelihood of finding H2 
is higher when AV is larger. 



7. Gas kinematics 

Figure[9]shows the distribution of the velocity widths of the low- 
ionisation metal line profiles. These measurements were com- 
pleted using the method described in Ledo ux et alj d2006al) . and 
many were previously published in that paper. It is clear from 
Fig. [9] that the probability of finding FL is higher when the ve- 
locity width of the low-ionisation metal line profiles is larger. 
Using the double-sided Kolmogorov-Smirnov test, we calculate 
a probability Pks = 0.075 that samples S and Shf are derived 
from the s ame parent distribut ion. A natural explanation is pro- 
vided bv lLedoux et al.l d2006al) who suggested that higher metal- 
licity DLAs arise from more massive objects. The amount of 
molecular gas and the st ar-formation rate in these systems could 
therefore be enhanced dHirashita & Fer rara 2005) naturally ex- 
plaining the above correlation. Another explanation could be 
that outflows can pull out cold and dusty gas, as well as provid- 
ing a large number of velocity components that wou ld increase 
the probability of finding H2 (Murp hy et al.l 120071) . However, 
H2 is usually detected in few components (typically one or 
two) in the main clumps of the systems, and not in the satel- 
lite components at high velocities that make up most of the 
width of the profile. The only exceptions to this are the DLA 



8. Evolution with redshift 

In Fig. [T0l we analyse the redshift distribution of DLAs within 
the UVES sample. The FL-bearing DLA sub-sample is sta- 
tistically indistinguishable from the overall UVES sample. A 
double-sided Kolmogorov-Smirnov test shows that sample S 
and sub-sample Sm have a probability Pks = 0.94 to be 
drawn from the same parent population. There is no indication 
of any evolution of the fraction of H2 -bearing DLAs with red- 
shift. Molecular hydrogen is more difficult to detect at the low- 
est redshifts observable using UVES (Zabs ^ 2) because fewer 
and weaker transitions are then covered by the UVES spectra 
and lines are located in the bluest part of the spectra where 
the signal-to-noise ratios are lower. While this could imply a 
bias in the detection rate of H2 in systems with low molecu- 
lar fractions, high molecular fractions (i.e., log/ > -4.5) are 
detected anyway because the corresponding column densities 
are always well above the detection limit. When considering 
only high molecular fractions, the same Kolmogorov-Smirnov 
test yields P ks(«SAShf) _ 0- 98. Contrary to previous tentative 
evidence dCurran et al.l 12004). there is also no evolution of the 
molecular fraction in systems with detected H2 (see Fig.lTOl). 
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However, because of small number statistics for systems at 
Z > 3, it is still unclear whether there is no really evolution in the 
sense of fewer H2 detections at higher redshift. Therefore, we 
restrict our claim of no evolution with redshift of the detection 
fraction and the molecular fraction to the range 1.8 < z a bs ^ 3. 
On the other hand, large molecular fractions are observed in the 
local Universe and significant evolutions of both the fraction of 
H2-detected systems and the molecular fractions are expected at 
lower redshifts. It is then of prime importance to cover the red- 
shift interval < z a b s < 1.8. This will soon be possible using the 
Cosmic Origins Spectrograph to be installed onboard the Hubble 
Space Telescope. 
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Fig. 10. Top panel: Redshift distributions of DLAs in the over- 
all UVES sample (S, 77 systems, solid histogram), systems with 
detected molecular hydrogen (<Sh2> 13 systems, grey), and sys- 
tems with log / > -4.5 (Shf, ten systems, dark grey). The inset 
shows the cumulative distributions from sample S (solid) and 
sub-sample Shf (dashed). The two distributions are statistically 
indistinguishable (Pks = 0.98). Bottom panel: Logarithm of the 
molecular fraction, log/, versus absorber redshift. There is a 
slight trend for more constraining upper limits at higher redshift. 
This is not physical however but due to the higher signal-to-noise 
ratios achieved towards the red, together with the larger number 
of H2 transitions covered by the spectra. 



9. Conclusion 

We present results of the largest survey of molecular hydrogen 
in high-redshift (1.8 < z a bs ^ 4.2) DLAs compiled to date 
using high signal-to-noise ratio, high spectral-resolution VLT- 



UVES data. We analyse data for 77 DLAs/strong sub-DLAs 
with N(Hi) > 10 20 cm" 2 , a dataset more than t wice as large as 
that studied previously by iLedoux et alj d2003l) . From the thir- 
teen high-redshift H2-bearing DLAs known to date, nine have 
been discovered by our group. Due to the superb quality of the 
Ultraviolet and Visual Echelle Spectrograph, we are able to de- 
tect unambiguously the H2 absorption features, measure accu- 
rate column densities, and in the cases of non-detections derive 
stringent upper limits. 

A double-sided Kolmogorov-Smirnov test shows that the ten 
H2-bearing systems with log/ > -4.5 (which is our conser- 
vative completeness limit) have Hi column densities that are 
compatible with those of the overall DLA population. This may 
be due however to small number statistics. There is evidence 
(see Sect. |4]l that the probability of finding large molecular frac- 
tions is higher in DLAs with large N(H 1), as observed in the 
LMC. About 7% of the systems with logjV(Hi) < 20.8 have 
log/ > -4.5 while ~ 19% of the systems with larger Hi col- 
umn densities have similar molecular fractions. There is no sharp 
transition in the molecular fraction of DLAs at any of the mea- 
sured total hydrogen (H 1+H2) column density. 

It is surprising to see that most of the DLAs have very low 
molecular fractions, i.e., log / ;S -6. This is much smaller than 
obser ved along li nes of sight in t he Galactic disk (ISavage et al.l 
1 19771) . However, W akke rj (2006) measured similar molecular 
fractions towards high-latitude Galactic lines of sight. 

We confirm that a good criterion to find H2-bearing DLAs is 
to select high-metallicity systems. Indeed, 35% of the systems 
with [X/H] > -1.3 have log / > -4.5 whilst this is the case for 
only 4% of those with metallicities lower than that. Since there is 
a correlation between metallicity and depletion factor, the latter 
being defined as [X/Fe] with X=Zn, S or Si, H2 is found in DLAs 
also having the highest depletion factors. Therefore, clouds with 
large molecular fractions a re expected to be dusty and clumpy 
(e.g., Hirashita et al. 2003). They can be missed because of the 
associated extinction and/or because of their small cross-section. 

The presence of H2 is closely related to the dust column den- 
sity. Indeed, all detections pertain to systems where the column 
density of iron into dust grains is larger than 5 x 10 14 cm" 2 , and 
about 40% of these systems have detectable H2. This shows that 
the presence of dust is an important ingredient in the formation 
of H2 in DLAs. The low molecular fractions measured for most 
of the DLAs are probably a consequence of low abundances of 
metals and dust. 

We show that the probability of finding H2 increases with 
increasing velocity width of low-ionisation metal line pro- 
files. The correlatio n between velocity width and metallicity 
(ILedoux et al.ll2006al) . if interpreted as a mass-metallicity rela- 
tion, provides a natural explanation in which H2-bearing DLA 
systems are preferably associated with massive objects where 
star-formation is enhanced. There is no evidence of systematic 
outflows in H2-bearing DLAs neither from the H2 profiles nor 
from the study of the associated high-ionisation phase. 

From the comparison between H2-bearing systems and the 
overall UVES sample, we show that there is no evolution with 
redshift of the fraction of H2-bearing DLAs nor of the molecular 
fraction in systems with detected H2 over the range 1 .8 < z a b s < 
3. This, compared to the large amounts of H2 observed in the 
local Universe, suggests that a significant increase of the molec- 
ular fraction in DLAs could take place at redshifts z a b s < 1.8. 
Ultraviolet observations from space are therefore needed to ob- 
serve the H2 Lyman and Werner bands in low and intermediate 
redshift DLAs. 



P. Noterdaeme et al.: Molecular hydrogen in high-z DLAs: The VLT/UVES database 



9 



Increasing the sample of H2-bearing DLAs is required to 
cover a large range in column density and derive the H2 fre- 
quency distribution, f(N(H2),X) - where X is the absorp- 
tion distance - in a way similar to studies completed for Hi 
dLanzetta et al.l 1 1 99 lb . This would assess whether the steep 
slope observed in f(N (Hi),X) at large column densities (e.g., 
Prochaska et al . 2005) is due to the conversion H 1— >H2 {Schaye 
2001) rather than to a magnitude-dependent bias from dust ob- 
scuration of the background quasars (e.g. , iBoisse et"ai] [l998; 
IVladilo & Perouxl20 05: Smett e et al.l2005l) . Both effects can ex- 
plain the apparent lack of DLAs of both high metallicity and 
large H 1 column density, as observed in all DLA samples. 

While the exact nature of DLAs is still an open debate, H2- 
bearing DLAs provide interesting probes of the physical condi- 
tions close to star-forming regions. About a decade ago, only one 
H2-bearing DLA was known. With the criteria revealed by the 
present survey, it will be possible to select these systems more 
efficiently and their numbers should increase rapidly. This opens 
up the exciting prospect of the detailed study of the ISM in dis- 
tant galaxies. 
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Table 1. Molecular and metal contents of UVES DLAs/sub-DLAs at 1.8 < z abs < 4.2. 
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< 13.8 


< 13.9 


< 


-6.15 


i 


Q2311-3721 


2.48 


2.182 


20.55 ± 0.07 


-1.57 ±0.10 


0.26 ± 0.05 


Si 


77 


< 14.7 


< 15.0 


< 


-5.00 


hj 


Q 2318-1107 


2.96 


1.989 


20.68 ± 0.05 


-0.85 ± 0.06 


0.42 ± 0.03 


Zn 


207 


15.49 


+0.03 
-0.03 

< 13.5 


a sa+o.08 
+-°"-nns 


f,g 


Q 2332-0924 


3.32 


3.057 


20.50 ± 0.07 


-1.33 ±0.08 


0.30 ± 0.03 


S 


111 


< 13.2 


< 


-6.45 


c 



Table 1. continued. 
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Quasar 
B1950 




Zabs 


log N(H I) 


[X/H] 


[X/Fe] 


X 


AV 
[kms- 1 ] 


logW(H 2 ) 
J=0 J= 1 


log/ 


Refs. 


Q 2343+1232 
Q 2344+ 1229 
Q 2348-0108 
Q 2348-0108 
Q 2348- 1444 


2.51 
2.76 
3.01 
3.01 
2.94 


2.431 
2.538 
2.426 
2.615 
2.279 


20.40 ± 0.07 
20.50 ±0.10 
20.50 ±0.10 
21.30 ±0.08 
20.63 ± 0.05 


-0.87 ±0.10 
-1.81 ±0.10 
-0.62 ±0.10 
-2.02 ± 0.08 
-2.10 ±0.08 


0.51 ±0.07 
0.18 ±0.01 
0.55 ± 0.02 
0.13 ±0.01 
0.25 ±0.10 


Zn 
Si 

s 

Si 

s 


289 
69 
248 
95 
55 


13 69 +u (w 

"■"'-0.09 

< 13.4 < 13.4 

18.52!}$ 

< 13.1 < 13.6 

< 13.6 < 13.6 


-6 41+ - 16 

-0.16 

< -6.40 

-1 69+ 038 

< -7.20 

< -6.38 


f,i 
h 

z, u 
h. u 
h 



" iLevshakov et al.l {2000, 200 H), b IPetitiean et al.l d2000h . c iLedoux et all d2003h, d iGe & BechtoTdl dl997T) , ' Petitiean eTall j2002h f 
iNoterdaeme et al.l d2007al), s [Sm ette et al. d2005t) , * this work, ' IPetitiean et al. J2006h . > lAkerman et alj d2005h, k ILevshakov et al.l d2002l) . ' 
ISrianand et all d2005t) , m iHeinmuller et al.l fcOOd), " ILevshakov & Varshalo 1 ^!] < 19851) , ° iFoltz et al] Jl988h, p Isrianand & PetitiearJ dl998h . 9 
ILedoux et all d2002h, r ISria nand et al. ( 200fj), iGe et alj d200lh . ' ILedoux et all d2006bh . " INoterdaeme et alj d2007bl) . 

+ A?(Fe n) from lProchaska et alj J2007bl). 

* Metal column densities from Prochaska et al. (2007b). Only an upper limit on A'(Feii) could be measured (Prochaska et al. 2003). The abundance 
of iron is therefore estimated through that of nickel, assuming [Fe/H] = [Ni/H] - 0. 1 dex. 



